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ABSTRACT 

We present the results of ^^CO J = 2 ^ 1 observations of the X-ray bright 
giant shell complex 30 Doradus in the Large Magellanic Cloud (LMC) using the 
Antarctic Submillimeter Telescope and Remote Observatory (AST/RO). This is 
the one of the largest HII complexes in the Local Group. We compare the ^^CO 
J = 2 ^ 1 observations against previously taken ^^CO J = 4 ^ 3 observations 
and analyze the spatial distribution of young stellar objects (YSO's) within the 
cloud using the Spitzer IRAC observations of the 30 Doradus complex. Both 
peaks of ^^CO J = 2 — > 1 and J = 4 — 3 emitting clouds coincide with the 
densest region of the filaments where multiple shells are colliding. We find that 
the YSO's are clustered in the southern ridge of the warm and dense molecular 
gas clouds traced by ^^CO J = 4 — > 3, indicating a filamentary structure of star 
formation throughout the 30 Doradus. We also find that some of Class I YSO's 
candidates which are likely to be Associated with a high-velocity component of 
^^CO J = 4 ^ 3 emitting clouds are present. This is a bona fide place where 
the triggered star formation had happened and newly formed stars may have 
produced such a high- velocity outflow interacting with the surrounding molecular 
cloud material. 
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1. Introduction 

The Large Magellanic Cloud (LMC) provides an excellent opportunity to study the 
effects of radiation from young star objects (YSO's) on their multi-phase interstellar medium 
(ISM). The LMC has a 27-degree inchne to the hne-of-sight and has small foreground and 
internal extinctions, making it possible to map the gas and dust of the ISM in the LMC 
without confusion with extraneous matter along the line-of-sight. The LMC stars are at a 
common distance from the Earth, unlike stars inside of the Milky Way's galactic plane, but 
are still close enough that individual stars and their stellar cjccta can be studied in great 
detail. We have chosen the 30 Doradus nebula inside of the LMC in particular given that Kim 
et al. (2005) detected a high- velocity molecular emission associated with the HII and HI 
shells, indicating that this is a very active region of star formation. Additionally, the LMC 
provides an excellent opportunity to study the effects of different UV radiation from stars on 
their environments in the multi-phase ISM, as well as to apply their observed mechanisms 
to studies of the early evolution of high-redshift, metal-poor galaxies. CO is the chosen 
species to trace the general distribution of molecular gas because it is the most abundant 
observable tracer of H2 and is sensitive to a density regime prevalent in the diffuse molecular 
gas surrounding dark clouds, molecular cores, and star-forming regions. For Galactic work, 
the ^^CO line is usually saturated and an isotope replacement molecule such as ^^CO is 
needed to obtain better information on column densities and cloud sizes. In the LMC, and 
therefore in 30 Doradus as well, the diffuse CO component is weaker due to lower abundance 
and higher photodissociation rates, so ^^CO is a less-biased molecular tracer which leads 
an assumption that the cloud is virialized and a reliable estimate of X conversion factor 
(Dickman 1978). 

The present study aimed to characterize the relationship between the warm and dense 
core of the molecular clouds traced with ^^CO J = 4 — > 3 and the distribution of young stars 
across the clouds in the 30 Doradus nebula. The mid- J CO emissions in the 30 Doradus 
region have recently been observed with the Antarctic Submillimctcr Telescope and Remote 
Observatory (AST/RO) and ^^CO J = 2 ^ 1 emission from the 30 Doradus reveal the 
extended structure beyond the ^^CO J = 4 3 emitting region. With the advent of 
the Spitzer Space Telescope, the distribution of young stellar objects across the region can 
be mapped with the Spitzer/IRAC photometry (Meixner et al. 2006). Jones et al. (2005) 
provide comprehensive study of protostars in the giant H II region complex N159 in the LMC 
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with the Spitzer IRAC observations. In this paper, color-color and color-magnitude diagram 
for Spitzer/IRAC sources have been used to study the nature of distributed populations and 
their locations in the molecular clouds in the 30 Doradus in conjunction with the mid- J CO 
emissions observed with the AST/RO. 

2. Observations 

The observations of ^^CO J = 2 — > 1 emission hues were performed during the aus- 
tral winter season of 2002 and 2003 at the Antarctic Submillimeter Telescope and Remote 
Observatory (AST/RO), a 1.7-m diameter, offset Gregorian telescope. AST/RO is located 
at 2847 m altitude at the Amundsen-Scott South Pole Station and capable of observing at 
wavelengths between 200 yum and 1.3 mm (Stark ct al. 2001). This site is very dry, so it is 
therefore a very good position from which to make submillimeter observations (Chamberlin, 
Lane, & Stark 1997). Since the observations of ^^CO J = 4 — > 3 emission line are described 
in Kim et al. (2005), we mention only the observations of ^^CO J = 2 — > 1 emission taken 
with the AST/RO in this paper. 

Emission from the ^^CO J = 2 ^ 1 line was mapped over a 30' x 30' region centered 
on R.A.=5'^38°^40% Dec.=-69°06'28" (J2000) with 90" spacing and a beam size of about 
180". The position switching mode was used and observation time at each position was 
typically 7 minutes and 20 seconds. The receiver used was a 230 GHz superconductor- 
insulator-superconductor (SIS) waveguide receiver with 79—90 K double-sideband (DSB) 
noise temperature. Two acousto-optical spectrometers (Schieder, Tolls, and Winnewisser 
1989) were used as backends. AST/RO had four available acousto-optical spectrometers 
(AOSs). An array AOS having four low-resolution spectrometer (LRS) channels with a 
bandwidth of 1 GHz (bandpass 1.6—2.6 GHz) and a resolution of 1 MHz corresponding to 
a velocity resolution of 1.3 km/s has been used for 230 GHz observation. LRS has 2048 
channels. AST/RO suffers pointing errors of the order of 1', and the beam sizes (FWHM) 
were 170-190" at 230 GHz. Atmosphere-corrected system temperatures ranged from 309 to 
364 K at 230 GHz. 

The standard chopper wheel calibration technique was employed, implemented at AST/RO 
by way of regular (every few minutes) observations of the sky and two blackbody loads of 
known temperature (Stark et al. 2001). Atmospheric transmission was monitored by regular 
skydips, and known, bright sources were observed every few hours to further check cali- 
bration and pointing. At periodic intervals and after tuning, the receivers were manually 
calibrated against a liquid-nitrogcn-temperaturc load and the two blackbody loads at ambi- 
ent temperature and about 100 K. The latter process also corrects for the dark current of 
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the AOS optical CCDs. The intensity calibration errors became as large as ±15% during 
poor weather periods. 

Once taken, the data in this survey were reduced using the COMB data reduction 
package. After elimination of scans deemed faulty for various instrumental or weather- 
related reasons (~ 10% of the total dataset), the data were Fourier transformed and the 
second highest amplitude channel was zeroed out, so well-defined baseline component was 
removed. Then, linear baselines were removed from the spectra. This is shown in Figure 1. 
The Fourier transforms are done using fast Fourier transform techniques and the spectrum 
is padded with zeros on the right to the next biggest power of 2 (e.g. 2048 channels). After 
a Fourier transform of the plotted channels of the original full spectrum were done and then 
selected high amplitude channel was removed, the data were re-transformed and inserted 
back into the original spectrum. Since the operations are all done in double precision, there 
is no loss of accuracy. 

In order to detect young stellar objects (YSOs), we used the Spitzer archival data of 
30 Doradus nebula (AORKEY 4379904) observed with the Infrared Array Camera (IRAC). 
IRAC data consist of four broadband images, each of which has a central wavelength of 3.6, 
4.5, 5.8, and 8.0//m (Fazio et al. 2004). First, we performed an array-location-dependent 
photometric correction on the Basic Calibrated Data (BCD) processed by the IRAC pipeline 
(version S13.2.0). This correction is required because the pipeline flat-fielding based on 
the zodiacal background is not appropriate for point source analysis. In the next step, 
we corrected "muxbleed" and "column pulldown" effects, which are image artifacts that 
appear around bright sources, and matched background levels of overlapping frames. Finally, 
mosaicked images were made from these corrected images using MOPEX mosaicker which 
is provided by the Spitzer Science Center (SSC). 

Source extraction and aperture photometry were performed using the Spitzer Astro- 
nomical Point Source Extractor(APEX) with 5 pixel aperture radius and 5-10 pixel back- 
ground annulus. Pixel phase corrections, which correct the dependence of photometry on 
the location of a source within its peak pixel, were applied on the photometry of channel 
1, and aperture corrections were applied on all channels. Then, we converted flux densi- 
ties into magnitudes using the IRAC zero-magnitude flux densities. All the values used 
for corrections and calculations were obtained from the IRAC Data Handbook version 3.0 
(http://ssc.spitzer.caltech.edu/irac/dh/ ). 
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3. Results 

We have observed ^^CO J = 2 ^ 1 line emission associated with the 30 Doradus nebula 
in the LMC (Figure 2). The peak of ^^CO J = 2 ^ 1 emission is only 40" away from the 
peak of ^^CO J = 4 ^ 3 emission. Including the pointing error and the beam size of ^^CO 
J = 2 — > 1 emission, the peak in ^^CO J = 4 — >^ 3 emission is very close to that in ^^CO 
J — 2 ^1 emission. Both peaks are not associated with the R136 cluster. Instead they are 
surrounded by RMS (WN6h), R139 (O6Iaf/WN)/R140 (WN6), and R144 (WN6h) multiple 
system north of the R136 cluster (see Figure 2). 

The peaks of ^^CO J = 4 — >^ 3 and ^^CO J = 2 — > 1 emissions are spatially distanced 
from the peaks of Ha and X-ray emissions as seen in Figure 2 and Figure 3 as it is expected. 
The apparent H II regions and the bright X-ray counterparts in this region indicate that both 
the fully ionized regions close to the massive stars, with the related ionizing radiation escaping 
from the HII regions, heat the H molecules and keep them dissociated. The morphology 
of Ha emission is almost identical to that of 8 /xm emission (Figure 4), tracer of Polycyclic 
Aromatic Hydrocarbon (PAH) emission, indicating that the PAH is only partially destroyed 
in the HII region. This is in contrast to the pretty good correlation shown in the SAGE 
study by Meixner et al. (2006) and SINGS study by Regan et al. (2006). High- velocity 
^^GO J = 4 — >^ 3 emitting gas detected in this region (Kim et al. 2005; Kim 2006) and the 
morphology of ^^GO J = 4 — >^ 3 emission indicate that the bulk of the GO emission must 
come from the entrainment of ambient atomic and molecular cloud gas by the stellar winds 
colliding with one another seen in the Ha image (Figure 2). Using the best estimate of the 
luminosity mass for ^^CO J = 4 ^ 3 outflow (Kim et al. 2005), 2.6±1.6x 10^ M©, and the 
dynamical age of the shells, 10 Myr, we estimate the rate of ambient gas entrainment to be 
Mfioyj=2.6 X lO^'^ Mq yr^^ with the product of preshock density and velocity, approximately 
nv ~ 18 cm~^ km s~^. The momentum injection rate is Pfiow = 1-3 x 10~^Mq km s~^yr~^ 
and the mechanical luminosity is to be Efioyj—5A x 10^ Lq. 

We detected the near-IR sources associated with the molecular cloud core in the ^^GO 
J = 2 — 1 and ^^GO J = 4 — > 3 emitting clouds in the 30 Doradus nebula by utilizing 
Spitzer GTO IRAG observations. We identified YSO candidates among these sources from 
the color-color diagram of 3.6 /xm — 4.5 /xm versus 5.8 /xm — 8.0 //m (Figure 5a), and 
classified them into three groups based on the criterion suggested by Allen et al. (2004). 
Since the reddening vector for 30 Doradus is about 1.3 mag, it can only change the colors of 
[3.6 //m] — [4.5 /im]=0.01 mag and [5.8 /im] — [8.0 ^m]=0.02 mag (Rieke and Lebofsky 1985). 
This is too small to affect the classification. These classes show the evolutionary sequence 
from Glass I to Glass HI. Glass I objects are surrounded by infalling envelopes. Glass II 
objects have accretion disks, and Glass I/II objects are in an intermediate state between 
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Classes I and II (Kenyon & Hartmann 1995, Megeath et al. 2004). They are denoted by 
open circles (Class I), squares (Class II), and diamonds (Class I/Class II) respectively in 
Figure 5. Class III objects were excluded from the classification because they have the SEDs 
of stellar photospheres and cannot be distinguished from foreground or background stars. 
Since the brightness of 8 ^m emission can be affected by PAH emission arising from the 30 
Doradus nebula, we present the color-color diagram of 3.6 //m — 4.5 /im versus 4.5 //m — 
5.8 iim in Figure 6. It confirms that the classified objects in Figure 5a fall in a relatively 
well-defined region of the color— color diagram in Figure 6. The distinction is clear. The 
only problem is that the 30Dor-09 object identified as Class II in Figure 5a shows much 
redder [4.5 yum] — [5.8 //m] color ~ 1.7. This object was identified as M red supergiant star. 
Dor IR 10 by McGregor & Hyland (1981) using the JHK photometry and IR spectroscopy. 
They also found 11 other M red supergiants in the 30 Doradus, but we could not find 
any corresponding sources associated with our classified YSO's catalog. Previous study by 
Brandncr ct al. (2001) reports 20 Class I protostars and Herbig Ae/Be candidates near 
R136 cluster in the 30 Doradus complex using JHK color-color diagram. 8 of these sources 
could be matched with the YSO population traced by present IRAC data analysis. Figure 
5b shows the location of these objects on the color-magnitude diagram of 3.6 jim — 8.0 //m 
versus 8.0 jim using the above symbols. Class I objects show strong 3.6 yt/m — 8.0 //m color 
excess with color indices of approximately 4, while most of Class II objects show relatively 
weak 3.6 jim — 8.0 /xm color excess with indices from 1 to 2. 



4. Discussion 

An HI aperture synthesis mosaicked map of the LMC with a spatial resolution of 50", 
created by combining 1344 separate pointings of the Austraha Telescope Compact Array 
(ATCA), shows an overall clumpy HI distribution featuring holes, shells, loops, filaments, 
and bubbles (Kim et al. 2003; Kim et al. 1998; Kim et al. 1999). The HI supergiant shells 
occupy a large volume of the ISM and a large number of giant shells, contained within each 
supergiant, are colliding with one another. The many shells observed often overlap and are 
interacting with one another, such as in the vicinity of the 30 Doradus complex where very 
active star formation has been happening simultaneously in many different centers. 

A number of the smaller shells formed on the rims of supergiant shells and are found in 
regions of very active star formation in the LMC (Kim et al. 1999). Recent submillimeter 
observations of the LMC supershells and supcrbubbles in the LMC using the Antarctic 
Submillimeter Telescope and Remote Observatory (AST/RO) aimed to look for the place 
where the triggered star formation had happened and found the ^^CO J = 4 — > 3 emission 
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toward the rim of the expanding giant HII complex in the 30 Doradus (Kim et al. 2005). 

They detected a possible high- velocity molecular emission associated with the HII and HI 
shells in the 30 Doradus nebula. This observational fact has been understood as the result 
of self-propagating star formation, where gravitational instabilities in the swept-up material 
of the supergiant shell caused fragmentation and a new round of star formation. Newly 
formed stars may have produced such a high velocity molecular outflows interacting with 
the surrounding molecular cloud material. Outflows and jets are commonly seen in the 
ambient medium of young stellar objects. 

It is commonly accepted that stars form in molecular clouds by the gravitational collapse 
of dense gas. Both low and high mass young stellar objects are embedded in the parent 
molecular and cause the heating of the surrounding gas. Often outflows from young stellar 
objects accelerate gas and are observed as molecular outflows. The radiative heating and 
magnetic fields are likely to be the main energy sources driving outflows for both low and 
high mass young stellar objects. Using the recent Spitzcr IRAC observations of 30 Doradus 
complex, we examine young stellar objects traced by IRAC data which might be associated 
with the high-velocity components of ^^CO J = 4 —> 3 emitting clouds. 

A total of 41 YSO candidates, which fall close to the core of the ^^CO J — 2 ^ 1 and 
^^CO J = 4 ^ 3 emitting clouds, were detected and classified from the IRAC survey of 
30 Doradus . We present our photometry results and classes in Table 1. Figure 4 presents 
the spatial distribution of detected YSO candidates towards the relatively dense molecular 
cloud core found in the ^^CO J = 4 ^ 3 and ^^CO J = 2 ^ 1 emitting gas. A strong 
association is not found between the location of YSO candidates and of dense molecular cloud 
core as seen in Figure 4. However, Class I candidates which are thought to be protostars 
surrounded by dusty infalling envelopes and exhibit signiflcant near-infrared extinction from 
their envelopes are very close to the peaks of ^^CO J = 4 — > 3 and ^^CO J — 2 ^ 1 emissions. 
Conspicuously, Class I YSO 30 Dor-30 is detected near the high-velocity component of ^^CO 
J = 4 ^ 3 emission. 13 YSO candidates are concentrated in the southern region of the ^^CO 
J = 2 ^ 1 emitting cloud, about 30 pc away from the molecular cloud core. Association 
densities are likely to be as low as 10^"^ stars pc^"^ at the distance of the LMC (Feast 1991) 
and correspond to the number density of molecular hydrogen n{H2) of 10^ molecules cm~^ 
(Kim 2006). It is notable that the relationship between the number of YSO candidates And 
their radial distribution from the core (Figure 4) reveals that more recently formed stars 
(Class I objects) are distributed more closely to the molecular cloud core. 

Figure 7 presents an analysis that measures the distance of the YSO's to the core of the 
CO cloud. The method performed was similar to one present in the paper by Teixeira et al. 
(2006), except that the measured variable is different. We utilized a Monte Carlo simulation 
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to analyze the difference between the observed spatial distribution and a randomly generated 
star field. We created 10,000 star fields with an equal number of randomly and uniformly 
placed YSO's in an equal observing area. The average distance to the cloud core of these 
objects was calculated and presented in a histogram to be compared with the distances of the 
observed YSO's. These comparisons show that there is a structural arrangement of YSO's 
with respect to the core of the CO cloud in the 30 Doradus Nebula, likely to be a filamentary 
distribution. In the total distribution of all classes, we find that there is a probabihty of 
0.15% of having more than 9 objects placed between 200 and 300 arcseconds from the core. 
This seems to be mostly due to the Class I YSO's given that there is a 0.02% probability 
of finding more than six Class I YSO's at that same distance. We can see that Class I/II 
YSO's are absent at this distance, but using the next bin (300-400 arcsec) we find that at its 
peak the Class 1/11 YSO's have a probability of 3.41% of having more than 4 objects. Class 
11 objects also follow the general distribution of the other two previous classes, but are too 
few and therefore tests are inconclusive. 



5. Conclusion 

AST/RO observations have revealed that the peaks of ^^CO J = 2 ^ 1 and ^^CO 
J = 4 ^ 3 emission are located at the rims of colliding H 11 shells (Chu and Kennicutt 1994) 
in the very active star-forming region in the 30 Doradus complex in the LMC. Previous 
studies by Kim et al. (2005) show that ^^CO J = 4 ^ 3 emission has a high-velocity 
component corresponding to the peak of the emission. This region must be the place where 
the triggered star formation had happened since newly formed stars may have produced such 
a high- velocity molecular outflow interacting with the surrounding molecular cloud material. 
Using the recent Spitzer IRAC observations of 30 Doradus complex, we find that Class I 
YSO's candidates are prone to be associated with the high-velocity component of ^^CO 
J = 4 ^ 3 emitting clouds. We also discover that the YSO's are clustered in the southern 
ridge of the relatively warm and dense molecular cloud cores traced by ^^CO J = 4 — 3, 
indicating a filamentary structure of star formation. 
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Table 1. YSO candidates near the 230 GHz ^^CO J = 2 — > 1 emission core in 30 Doradus 

complex 
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*These sources correspond to the YSO candidates studied by Brandner et al. (2001). 




100 150 200 250 300 350 

V^, (kms-0 

Fig. 1. — The ^^CO J = 2 1 raw spectra (Top: integration of 1 min and 30 sec), Fourier 
Transformed and baseline removed spectra (Middle), and integrated for 7 min and 20 sec 
spectra (Bottom) observed towards the peak of ^^CO J = 4 — >■ 3 emission in the 30 Doradus 
complex. 



- 13 - 




Right Ascension (J2000) Right Ascension (J2000) 



Fig. 2. — The ^^CO J = 4 — > 3 emission line (blue contour) is displayed and compared to 
the ^^CO J = 2 1 emission line (red contour) from this position. These are overlayed on 
a grayscale image of Ha. The contour levels are 1.3, 1.7, 2.1, 2.5, 2.9, and 3.3 K km/s. Ha 
image is taken with a CCD camera mounted at the Curtis Schmidt telescope at CTIO by C. 
Smith. 



-14- 




• 



s'^ag^ao^ 39^00'' 38™30'' 

Right Ascension (J2000) 

Fig. 3. — The ^^CO J = 4 ^ 3 and ^^COJ = 1 — > emission lines are overlayed on the 
Chandra X-ray image of 30 Doradus. Chandra ACIS-I 0.5-10 keV observations were exposure 
corrected and the resultant image was obtained from binning by 4 pixels and has angular 
resolution of 2" . 
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Right Ascension (J2000) 

Fig. 4. — Molecular cloud core found in the ^^CO J = 2 1 emission map is overlaid on 
the 8.0/im Spitzer image. The 41 YSO candidates are marked by the same symbols used in 
Figure 5. 
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Fig. 5. — (a) Color-color diagram of 3.6 fim - 4.5 iim versus 5.8 fim - 8.0 fim. Detected 
infrared sources are classified based on the criterion suggested by Allen et al. (2004) and 
Megeath et al. (2004). Open circles represent Class 1 objects and diamonds represent Class 
1/11 objects. Squares represent Class 11 objects and Filled circles represent both Class 111 
objects and ordinary stars. Solid box represents the region of Class 11 sources (Allen et al. 
2004). The solid hne shows the division between the Class I and Class I/II candidates, (b) 
Color-magnitude diagram of 3.6 /im - 8.0 /im versus 8.0 /im using the same symbols. 
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Fig. 6. — Color-Color Diagram of 3.6 //m - 4.5 //m versus 4.5 yum - 5.8 yum. The 41 YSO 
candidates are marked by the same symbols used in Figure 5. We plot the locations of the 
late spectral types from Jones et al. (2005) along the vertical axis. 
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Class II YSO's Detected in 30 Doradus 
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Fig. 7. — The numbers of detected objects for each class are plotted against the radial 
distance from the ^^CO cloud core. The sohd lines are the objects observed with the 
IRAC/Spitzer data, while the dotted hues represent the average values of the randomly 
generated fields. 



